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ABSTRACT 

We perform spherically- symmetric general-relativistic simulations of core collapse and the postbounce pre- 
explosion phase in 32 presupernova stellar models of solar metallicity with zero-age-main- sequence masses of 
12 Mq to 12OM0. Using energy-dependent three-species neutrino transport in the two-moment approximation 
with an analytic closure, we show that the emitted neutrino luminosities and spectra follow very systematic 
trends that are correlated with the compactness M/R) of the progenitor star's inner regions via the accretion 
rate in the preexplosion phase. We find that these qualitative trends depend only weakly on the nuclear equa- 
tion of state, but quantitative observational statements will require independent constraints on the equation of 
state and the rotation rate of the core as well as a more complete understanding of neutrino oscillations. We 
investigate the simulated response of water Cherenkov detectors to the electron antineutrino fluxes from our 
models and find that the large statistics of a galactic core collapse event may allow robust conclusions on the 
inner structure of the progenitor star. 

Subject headings: equation of state - hydrodynamics - neutrinos - stars: evolution - stars: neutron - stars: 
supemovae: general 



1. INTRODUCTION 

The radial instability of the electron-degenerate, 
Chandrasekhar-mass core marks the beginning of the fi- 
nal episode in the life of a massive star with zero-age 
main-sequence (ZAMS) mass in the range ^ SM© - ISOMq. 
Collapse ensues and, once fully dynamical, separates the core 
into a subsonically homologously contracting inner core and 
a supersonically collapsing outer core. At nuclear density, the 
repulsive component of the nuclear force leads to a stiffening 
of the equation of state (EOS). This stabilizes the inner core, 
which overshoots its new equilibrium, then rebounds into 
the outer core, launching a strong hydrodynamic shock wave 
from its edge. This instant in time is referred to as core 
bounce. The inner core has a mass of ~ 0.5 at bounce and 
this material becomes the unshocked core of the protoneutron 
star. The shock formed at core bounce propagates into the 
outer core, but dissociation of accreting nuclei into neutrons 
and protons and electron capture on free protons in the region 
behind the shock (the postshock region) soon sap its might, 
driving it into submission to the ram pressure of accretion. 
The shock stalls and turns into a standing accretion shock 
that must be revived to unbind the stellar envelope and drive 
a core-collapse supernova explosion. 

Neutrinos play a pivotal and dominant role in stellar col- 
lapse and core-collapse supernovae. Neutrinos and antineu- 
trinos of all flavors carry away the ~300B (= 3 x 10^^ ergs) 
of gravitational binding energy of the remnant neutron star 
over tens of seconds after core bounce. Aided by multi- 
dimensional fluid instabilities, they probably deposit, within 
a few hundred milliseconds after bounce, sufficient energy in 
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the region behind the stalled shock to revive the sho ck and 
make typical ^ 1 B core-colla pse supernova explosion ( [Janka 
'et al.'2007!'Mul ler et aL]j 20l2 and references therein). Only 
hyper-energetic (i.e ., (^(IQ) B) explosion s may require a dif- 
ferent mechanism ( [Ugliano et al. 2Q12| ), e.g., rapid rotation 



combined with strong magnetic fi elds, which may lead to en- 
ergetic jet-driven explosions (e.g., [Burrows et al.|2007| ). 



For a galactic or near-extragalactic core-collapse super- 
nova, neutrinos offer the unique possibility of directly ob- 
serving the dynamics and thermodynamic conditions preva- 
lent in the supernova core. Together with gravitational waves 
(see, e.g., Ott 2009; Kotake 2011) they wiU herald the next 
nearby supernova possibly hours before any telescope sensi- 
tive to electromagnetic waves will notice the event. In the 
probable case that the next galactic supernova occurs in a 
dust-enshrouded region and/or close or behind the galactic 
center, the supernova may be impossible to observe in broad 
bands of the electromagnetic spectrum, making neutrino and 
gravitational- wave observations even more important. 

The observation o f neutrinos fr om SN 1 9 87 A in the Large 
Magellanic Cloud ( Hirata et aL] p'987: Bi onta et aL][T987 



Alekseev et al. 1987| confirmed the basic picture of core 



collapse and earTy"protoneutron star evolution (e.g., Sato & 



'Suzuki"" 1987! "Bruenn"'1987'; 'Burrows & Lattime r||1986J|Bur- 



rows 1987, 1988; Arnett et al. 1989; Jegerlehner et al.|1996 



Loredo & Lamb 2002; Yuksel & Beacom 2007; PagfiaroTil 
et al. 2009b and references therein), but the small number and 
poor timing of the observed events did not allow far-reaching 
and robust conclusions on core-collapse supernova dynamics 
and the involved neutrino physics, nuclear physics and astro- 
physics. 

The situation will be completely different when the 
neutrino burst from a galactic core-collapse supernova 
reaches current and near- futu re neutrino dete c tors on Earth. 
Super-Kamiokande (Fukuda et al.| [2003 [ [Ikeda et al. 
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[Super -Kamiokande Collaboration] ' '2007 ), IceCu be (Ab- 
basi et al. [IceCube Collaboration] 20lTj) LVD ( [Aglietta 
eTaf. 1992yyigorito [LVD Collaborationf'MTr), Borexmo 
lAlimon ti et alT [B OREXINQ Collaboration] 2009; Cado- 
nati et al.||2002D KamLAN D (|Piepke||2001|), SNO+ (|Kraus 
& Peeters| |20TQ |), Noi^a (D avies, [for the NOz/A Collabo- 



ration]! [2011 1 ) and others will together see many thousands 



of neu trinos from a core collapse event at lOkpc ([Scholberg 



|2012|). Distance estimate s based on the obs erved neutrino flux 



(KachelrieB et al."2005), sky localization ( [Beacom & Vogel| 
P"999 i Tomas et al. 2003 ), and triggering of gravitational- wave 
searches by the reconstruct ion of the time of c ore bounce 
( [Pagharoli et"aLl[2009a, ,Halzen & Raffelt|[20Q9l ) will likely 
all be possible. 

Well-timed high-statistics coincident neutrino observations 
will allow to probe in detail a broad range of supern ova astro- 



physics, nuclear physics, and neutrino physic s (see [Burrows 
|etaL][T992l |Wurm et al.||2QT2l |Raffelt||2QTQl for overviews). 
Fast characteristic temporal variations in the preexplosion 
neutrino fluxes would be tell-tale signs of multi-dimensional 
fluid instabilities in the postshock reg i on (Qtt et al.,,2008[ 
IMarek et al.| [20091 |Lund et al.|[2QTQl [Brandt et al.[[2QTTt 
and/or early postbounce rin g-down oscillat ions of a rapidly 
spinning protoneutron star ( Ott et al.][2012| ). A sudden deep 
drop of the accretion-driven component of the neutrino lu- 
minosity (primarily in and Vq) within a few hun dred mil- 
liseconds would indicate the onset of explosion (e.g., [Burrows 



eraLl[T992l [Hudepohl et ^^[MTqI [Fischer et al.[|2012| [2010 
and an abrupt cut-off of the entire neutrino flux within sec- 
ondsof its onset would indicate the forma tion of a black hole 
jBuTOws 1986, 1988; Beacom etal. 2001 ;ILi ebendorfer et al.| 
|2003, 2004; Sumiy oshi et al . 2006; Fischer et al."2009). The 
spectral characteristics and long-term spectral evolution of the 
neutrino flux could provide important constrai nts on the nu- 
clear EOS ( ^Roberts et a l.'2012VMarek et al.'2009) and/or the 
spin of the progenitor core ( Ott et al. 2008; Marek & Jankaj 



2009| ). The ability of some detectors to distinguish interac- 



tions of different neutrino flavors would lead to constraints on 
the neutron-to-proton ratio in the neutrino-driven wind phase, 
allowing an observational test of core-collap se supernovae as 



potential sites for r-process nucleosynthesis ( jHudepohl et al. 
|20T0HFischer et al.|2010||Wurm et al.|2QT2l ). 

The neutrino signature of core collapse, of the subsequent 
core-collapse supernova evolution and of the protoneutron 
star cooling phase, is invariably intertwined with neutrino os- 
cillation physics. The robustness of all of the above men- 
tioned observational conclusions will depend on our under- 
standing of the impact of neutrino flavor oscillations. Neutri- 
nos propagating from their emission site to detectors on Earth 
may experience (/) so-cal led vac uum oscillations driven by 
neutrino mass differences ([ Pont ecorvo 1968 ), (//) oscillations 
mediated by a resonance in v-e~ scat tering (t he Mikheyev- 
Smirnov-Wolfenstein [M SW] effect; [Mikheev & Smimov] 
|1985t [Wolfenstein|[1978j), an d (/ //) oscillations du e to v-v 
scattering (Pantaleone 19921 see Duan et al. 2010 1 for a re- 
view). Vacuum and MSW oscillations are well understood 
and their outcomes depend essentially only on neutrino mix- 
ing parameters, in particular the neutrino mass hierarchy and 
the mixing angles. The z^-z^- scattering driven oscillations, on 
the other hand, have a non-linear Hamiltonian that may lead 
to so-called collective oscillations with very complex spatial 
and temporal outcome that remains to be fully understood 
(see, e.g. „Hannestad et al.|2006,puan et al.|2007[[Fogli et al.| 



20091 [Dasgupta & Dighe||20Q8l IDuan et al.[|2QT 0l and refer- 
ences therein). However, a number of recent studies suggest 
that collective oscillations may be completely or at least par- 
tially suppressed in the preexplo sion accretion phase of or- 
dinary core-collapse superno vae ([Chakr aborty et al.'2011b'a* 
Sarikas et al. 2012 ), but see Cherry et al.| ( ,2012j and ,Dasgupta| 
et al.[ ( [2012| ) for discrepant results. 

Provided that collective oscillations can be ignored in the 
preexplosion phase and that the ^13 mixing angle i ndeed has 
the la rge value suggested by recent measurements ( [An et al. 
2012| ), the neutrino mass hierarchy may be inferred from 
the qualitative shape of the early postbounc e neutrino signal 
( [KachelrieB et al.|2005t[Serpico et al.|2012[ ). 

A preexplosion accretion phase with suppressed collective 
oscillations would also offer the opportunity to probe the 
structure of the progenitor star on the basis of the observed 
neutrino signal. The details of the preexplosion neutrin o emis- 
sion h ave been discussed carefully, e .g., by Thomps on et al. 
( [20031 ) and [Liebendorfer et al.[ ( [2004j ) and we shall not repeat 
them here. It is, however, necessary to outline its most salient 
features. For simplicity, we neglect neutrino oscillations in 
the following. 

In core collapse and in the subsequent postbounce evolu- 
tion, emission of and occurs via charged and neutral 
currents, while heavy-lepton neutrinos - {z^/^, ^/x, ^r, ^r} 
are created exclusively via thermal neutral-current pair pro- 
cesses. Before core bounce, only are emitted from elec- 
tron capture in the collapsing core. Milliseconds after core 
bounce, the shock breaks out of the neutrino sphere (where 
the optical depth is r^^ ^2/3) and a strong burst of is 
emitted for '^20 ms from rapid electron capture on the freshly 
abundant free protons behind the shock, are copiously cre- 
ated in the hot interior of the protoneutron star after bounce 
and begin to diffuse out, leading to a steep rise, quick lev- 
eling and subsequent slow decay of the luminosity (L^J. 
Ve production via charged-current positron capture is initially 
suppressed due to the high degeneracy of the electrons. The 
latter is partially lifted after bounce at the moderate-density, 
hot edge of the protoneutron star and Ly^ rises, reaching or 
surpassing the value at which Ly^ levels off after the neutron- 
ization burst decays. The subsequent preexplosion luminosity 
can roughly be split into a diffusive component from the core 
and accretion luminosity (ex GM\Ry\M jR^, where is an 
approximate ne utrino sphere rad ius) from or from above the 
neutrinosphere jBurrows 1988). L^^ is primarily diffusive, 
while Ly^ and Ly^ are dominated by accretion. In general, 

l^V, ~ l^V, > Ly,, but 4Ly^ = Ly^ + Ly ^ + Ly^ +Ly^ > Ly^ + Ly^ . 

Vx have the lowest opacity, since they interact only via neutral 
currents. They decouple from matter at the smallest radii and 
highest temperatures and thus have the highest average ener- 
gies {ey). Ve have a slightly lower opacity than v^, leading to 
the well established neutrino energy hierarchy in the preex- 
plosion phase {ty^ > (cy^) > (e^J. The mean energy of all 
species grows with increasing postbounce time, reflecting the 
recession of the neutrinospheres due to the contraction of the 
protoneutron star. 

Considering that the accretion luminosity will scale with 
the postbounce accretion rate M, one would naturally expect 
an increase of the detected neutrino events with increasing 
mass of the stellar core. Since higher accretion rates corre- 
spond to more material compressing and settling more rapidly 
on the protoneutron star, the latter' s outer regions will be hot- 
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ter. Thus the thermal neutral-current emission will be en- 
hanced, leading to higher luminosities and higher mean neu- 
trino energies. 

The variation of the preexplosion neutrino s ignal with pro- 
genitor star ZAMS mass was first discussed by |Woosley et al. 
p986) based on the pioneering simulations of Wilson|^1985 



and Wilson et al. (1986). These authors provided total emis 
sion characteristics and spectra that show a systematic in- 
crease of total energy emitted in neutrinos and mean en- 
ergy w ith ZAMS mass in the range from 1O-25M0. Mayle 



|et al.| ( ^198 7 ), before SN 1987A, carried out simulations of 
a range of progenitor stars with ZAMS mass in the range 
12-1OOM0. They found that the neutronization burst 
shows little dependence on the progenitor, due to the rather 
universal homologous collapse and bounce dynamics. Fur- 
thermore, they mentioned, though did not discuss in detail, 
that the luminosities and mean neutrino energies increase as a 
function of iron core mass (and not ZAMS mass). A more de- 
tailed and clear physical discussion was provided by Bruenn 
|1987), who contrasted the predicted neutrino signal from 
the early postbounce phase in two different progenitor core 
models with neutrino observations of SN 1987 A. He noted 
that there are significant uncertainties in connecting a given 
ZAMS mass to precollapse structure. Instead of a progeni- 
tor with an associated ZAMS mass, he considered a massive 
(and high-entropy) 2.O5-M0 iron core model and a lower- 
mass (and lower-entropy) 1.35-M0 iron core model in his 
spherically-symmetric (ID) neutrino radiation-hydrodynamic 
simulations. He showed that the more massive core leads to a 
consistently higher luminosity in both the accretion and dif- 
fusion sectors. The water Cherenkov detectors that observed 
neutrinos from SN 1987 A are most sensitive to the in verse 
beta decay (IBD) reaction v^^-p^ ^ + Bruenn ( 1987 ) pre- 
dicted a factor of two difference in the integrated early IBD 
events between the massive and the low-mass core in these 
detectors. He concluded that the neutrino signal observed by 
these detectors from SN 1987A was most consistent with the 
low-mass core. |Burrows| (p^88 ), who carried out a parameter 
study of quasi-hydrostatic protoneutron star cooling, consid- 
ering various initial masses, ad-hoc accretion rates, and dif- 
ferent nuclear EOS, found a similar trend. He showed that 
more massive cores, higher accretion rates, and softer EOS 
lead to stronger, higher-energy neutrino emission. Some of 
his strongest emitters were cases in which eventually a black 
hole was formed. 

Liebendorfer and collaborators carried out a sequence of 
studies of the progenitor dependence of the neutrino signal us- 
ing modem general relativistic ID radiation-hydrody namics 
simulations ( [Liebendorfer et al.|[2QQT] [20021 [20031 [2004 ). 
They showed that the neutronization burst is indeed al- 
most independent of progenitor structure (as first suggested 
by Mayle et al. 1987 ). They also qualitatively and quantita- 
tively connected the evolution of the postbounce preexplosion 
luminosity to the postbounce accretion rate, but did not dis- 
cuss observational implications. Their results were corrobo- 
rated by similarly sophisticated subsequent studies of Kachel- 
rieB et al.| (|200 5); Bura s et al.| (|2006|); iFischer et a l. (2009); 
Sumiyoshi et al. (2 008); [Fischer e t al.li2 010[[2012 ); Serpico 
et al. (2012). Of these, [BuraseTa l. (2006) presented the most 
comprehensive analysis and also compared between ID and 
axisymmetric (2D) results. They found that in 2D, convection 
in the protoneutron star alters the structure of the latter, affect- 
ing the neutrino emission starting ~ 100 ms after bounce, but 



preserving th e overall systematics with accretion rate. Buras 
et al.[ p006| ) also were the only authors to suggest that the 
accretion-rate dependence of luminosity and total emitted en- 
ergy in the preexplosion phase could be used to infer the struc- 
ture of the progenitor. The other studies, being focused on 
aspects such as neutrino oscillations, black hole formation, or 
the late-time post-explosion evolution, did not consider obser- 
vational consequences. 



[Thompson et al.[ ( [2003 1 ), using a limited set of three pro- 
genitor models ({11, 15, 2O}M0 at ZAMS), found similar 
systematics as the aforementioned studies, but also carried 
out an analysis of the expected signal in various neutrino de- 
tectors in the first 250 ms after bounce. They computed IBD 
event rates for their 2O-M0 and II-M0 and found a factor 
of two more IBD events for the former, which would allow 
a high-confidence distinction between these progenitors for 
a galactic core collapse event. However, their 15-M0 model 
yielded a postbounce neutrino signal very similar to that of 
their II-M0 model and would be indistinguishable by neu- 
trino observations alone. This suggests that ZAMS mass is 
not a go od parameter to describe presupernova stellar struc- 
ture (cf., [Bruenn|1987l ). 

In this article, we present a fresh look at the progenitor de- 
pendence of the neutrino signature in the preexplosion accre- 
tion phase of core-collapse supernovae. We perform ID gen- 
eral relativistic radiation-hydrodynamics core collapse sim- 
ulations of 32 progenitor models from t he single- star solar- 
metall icity presupernova model suite of [Woosley & Heger 
( [2007| ) and follow the postbounce preexplosion evolution for 
450 ms. In ZAMS mass, these models range from 12 Mq 
to I2OM0, but guided by the previous results discussed in 
the above, we choose not to parameterize our simulations by 
ZAMS mass. Instead we employ the compactness parameter 
^ M/R(M) (for a relevant mass scale M, measured a t the 
time of bounce). As shown in [O'Connor & Ott[ ( [201 1| ), .^m 
is the most important stellar structure parameter governing 
the postbounce accretion evolution to a remnant mass scale 
M. We demonstrate that the preexplosion neutrino emission is 
very well parameterized by the compactness. The preexplo- 
sion luminosities and mean energies of all neutrino species 
increase essentially monotonically with increasing ^m- We 
compute predicted integrated IBD events for a galactic core- 
collapse supernova in the Super-Kamiokande detector and 
show that the clear systematics governed by carries over 
to observation, even when standard MSW neutrino oscilla- 
tions are taken into account. Our results thus indicate that 
- in the absence of complicated collective neutrino oscilla- 
tions - a high-statistics detection of neutrinos from the pre- 
explosion phase will allow, in principle, a tight constraint of 
the compactness of the progenitor star's core. This, however, 
will require knowledge of the nuclear EOS and of the rotation 
rate of the collapsed core, since, as we show, both can dilute 
the otherwise clear compactness-dependent neutrino emission 
systematics. 

This paper is structured as follows. In ^ we discuss our 
general relativistic hydrodynamics code GRID and introduce 
its extension to neutrino radiation-hydrodynamics in the two- 
moment approximation, nuGRlD. The initial models and the 
employed EOS are discussed in ^ In ^ we present re- 
sults from a benchmark collapse and postbounce simulation 
that allows us to compare with the previously published code 
comparison of Liebendorfer et al.,(2005) to assess nuGRlD's 
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ability to reproduce results of full Boltzmann neutrino trans- 
port. We present the results of our simulations in ^ analyze 
the dependence of the neutrino signal on progenitor compact- 
ness, discuss predicted IBD signals from a galactic core col- 
lapse event in the Super- Kamiokande detector, and explore 
potential degeneracies introduced by EOS and rotation. Fi- 
nally, in ^ we critically summarize our work and conclude 
by contrasting our results with the early neutrino signal ob- 
served from SN 1987A. 

2. METHODS 

We make use of the open - source ID general rela tivistic 
hydrodynamics code GRID ( [O'Connor & Ott|[2QTm avail- 
able at http : / / www .stellar coll apse, org) outfit- 
ted with an energy-dependent multi-species Ml neutrino 
transport scheme in which the zeroth and first moments of 
the neutri no distribution function are evolved. We refer the 
reader to [O'Connor & Ott| ( [2QTQl ) for details on GRID and 
describe in the following our current implementation of the 
transport scheme. For this first application, we neglect the 
computationally-expensive energy-coupling neutrino interac- 
tions and transport terms - these terms are undoubtedly im- 
portant for making highly accurate predic tions of the neutrino 
signature (see, e.g., [Lentz et al.|2Q12a|b| ), but are unlikely to 
affect the general trends we observe. We will address them 
in future work, but provide a discussion on the consequences 
of neglecting these terms via a comparison to full Boltzmann 
neutrino transport simulations in ^ 

Our Ml scheme closely follows Shibata et al.jpOll] ), who 
formulate the Ml evolution equations in a closed covariant 
form. The scheme is simplified greatly by neglecting the 
energy-coupling terms. This further requires that the veloc- 
ity dependent terms are also ignored. In this limit, and us- 
ing the Schwarzschild-like metric and radial-polar slicing of 
GRID and setting G = c = Mq = 1, the evolution equations for 
the neutrino energy density, E(^y), and the neutrino flux vector. 



FrXu), s implify from Equations 3.37 and 3.38 of Shibata et al. 

( |2mT] )to 



a 



(1) 



and 



X2 



where 5" is the neutrino interaction source term (see below), 
a is the lapse function and X = (1 -2M(r)/r)~^/^. Prr,(u) is the 
neutrino pressure tensor and is taken to be an interpolation be- 
tween the two limiting cases of free streaming and diffusion. 
We follow Shibata et al. (2011 ), who express Pu^u) ^.s 



_ 3/?(.)-l , 3(1 -P(.)) 

~ 2 ^n,(z^),thin M/,(z^),thick ? ) 

where p(^j,) is the Eddington factor, taken here to be the max- 
imum entropy closure in a closed, a nalytic form ( [Minerbo 
T978l[Cernohorsky & Bludman|1994| , 



(4) 



In the no-velocity limit for GRID, /(j,) = |iV(j,)/(£^(j,)X)|. The 
free streaming and diffusion limits of the neutrino pressure 



tensor are PrrMM^ = E(j,)X^ and Pmu),Mck = giiE(u)/3, respec- 
tively. 

We set out to solve the system of equations via standard hy- 
perbolic methods borrowed from conservative hydrodynamic 
schemes ( [Pons et aL]|2QQ0[ ). Transport variables live at cell 
centers and we employ piece-wise l inear reconstruction to 
cell interfaces with van Leer's limiter ([van Leer 1977 ) and the 
HLLE approximate Riemann solver (Einfeldt 1988 ) for calcu- 
lating the intercell fluxes. A complication arises when solv- 
ing the neutrino moment equations in the high-opacity limit. 
In this case, the standard fluxes returned from the Riemann 
solv er are domin ated by a numerical diffusion term. We fol- 
low [Audit et aL[ (|20Q2) and modify the fluxes to correct for 
this. Essentially, tnemodified fluxes return the diffusion-limit 
flux in the high opacity limit. The modified form of the neu- 
trino energy density flux through the / + 1/2 interface is given 
by 



(5) 



The corresponding modified interface flux for the neutrino 
flux evolution equation is. 



= ^^P^ + (1 - +^^))/2 , (6) 



with 



a^-a 



(7) 



In these equations, e(j,) controls the modific ation to the 



fluxes to account for the high opacity. Following [Audit et aL 
( |2002| ), we take 



e(^) = min 1, 



1 



(8) 



where hz^jy) is the sum of the scattering and absorptive opac- 
ities. These opacities are strong functions of energy and are 
also species dependent. We note that when e^^) is 1, the inter- 
cell fluxes reduce to the standard HLLE approximation. The 
characteristic speeds needed for the HLLE scheme are ca lcu- 
lated in the same spirit as the ne utrino pressure tensor ( Shibata[ 
[et al.|2Qllt[Kuroda et al.|20T2l ), 

_ 3/7(^)-l , 3(1 

\u) - 2 ^(z^),thin H ^(zy),thick , {y) 

where in the zero velocity limit, A(j,)^thin = and 
A(i/),thick = ±a/(>/3X). and a' are the maximum and mini- 
mum values, respectively, of these characteristic speeds eval- 
uated from both the right and left reconstructed variables. 

Finally, the source terms in Equations [l] and |2| are taken 
from jShibata et al.[ ( [2011[ ). In the zero velocity limit. 



(10) 
(11) 



where /^(^y), tia,{v)^ and ^s^y) are the neutrino emissivity, neu- 
trino absorption opacity, and the neutrino scattering opac- 
ity, respectively. We precompute the neutrino interaction 
terms for each neutrinos species (we treat and = 

{^/27^/i7^r7^r}) and neutrino energy group in dense tabular 
form as a function of density p, temperature T, and electron 
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fraction Yg. We then use linear interpolation for efficient on- 
the-fly interpolation. We include all standard iso-energetic 
scattering processes, charged-current absorption and emis- 
sion, and thermal p air-production processes (Burro ws et al. 
|2Q06t |Bruenn|[l985 ) in the calculation of the neutrino in- 
teraction terms. Since the neutrino-matter interactions for 
heavy-lepton neutrinos and antineutrinos are slightly differ- 
ent, NuLib averages the two values of the emissivities and 
opacities. Our library of neutrino interaction routines (which 
we call NuLib) is open source and available as a GitHub 
repository at http : / /www . nulib . org. NuLib requires 
an EOS for the evaluation of the emissivities and opacities. 
Our treatment of thermal pair processes in GRID warrants 
some comments. Since we do not currently consider energy 
(or species) coupling for thermal emission processes such as 
electron-positron annihilation to a neutrino-antineutrino pair, 
we compute an emissivity based on the thermal content of the 
matter ignoring any final state neutrino blocking. We limit the 
neutrino energy density to the blackbody occupation density 
by using Kirchhoff 's law to determine an effective opacity for 
neutrino-antineutrino annihilation from the thermal emissiv- 
ity. As we shall see, this method performs well at predicting 
the thermal neutrino flux of the heavy-lepton neutrinos during 
the preexplosion phase. 

In nu GRID, we first update the hydrodynamic variables to 
the ^+ 1-th timestep. We then compute the neutrino opacities 
and emissivities associated with the updated hydrodynamic 
variables. We update the radiation field operator- split. The 
flux term is solved explicitly, using the radiation moments of 
the ^-th timestep. We calculate the neutrino-matter interac- 
tion terms using the n-\-l radiation moments via a local im- 
plicit update. With the ^ + 1 radiation energy density source 
term, we then update the energy density and electron fraction 
of the matter. We use 24 energy groups, with lowest-energy 
group centers at 0.5 MeV and 1.5 MeV, and then spaced loga- 
rithmically up to 200 MeV for Ue, Ve and v^. We note that for 
the highest energy bins it occasionally occurs that the evolved 
neutrino flux vector exceeds the evolved neutrino energy den- 
sity. This tends to occur in the most dynamic phases of our 
simulations and where the opacities vary significantly from 
one zone to the next. When this is the case we limit the neu- 
trino flux to the neutrino energy density. 

3. INITIAL MODELS AND EQUATIONS OF STATE 

We employ the most recent non-rotating solar-metallicity 
single-star mo del set from the stellar evolution code KEPLER 
( [Woosley & H eger 2007). This model set contains the pre- 
supemova configuration of 32 stars ranging in ZAMS mass 
from 12M0 to 12OM0. We denote individual models by 
^XXWHOT, where XX corresponds to the integer ZAMS mass 
of the model, e.g., ^12WH07 is the 12-M0 model of this 
model set. In O'Connor & Ott| (2011), we investigated this 
and other model sets in the context of black hole formation. 
Under the assumption of a failed core-collapse supernova, we 
found a strong empirical relation between the properties of 
the presupernova structure and the evolution of the failing su- 
pernova, e.g., the time to black hole formation. This led to 
a clear prediction: If we observe black hole formation in a 
failed core-collapse supernova via neutrinos, the lifetime of 
the protoneutron star (and thus of the neutrino signal) relays 
direct information about the presupernova structure. How- 
ever, such a prediction, (/) requires a failed supernova, which 
may not be the norm, and (//) has a strong dependence on the 
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Fig. 1. — Compactness parameters for the 32 considered presupernova 
models of Woosley & Heger (2007) versus ZAMS mass as evaluated from 
collapse simulations with the LS220 EOS. We show both ^1.75 and ^2.5- 
The mapping between ZAMS mass and precollapse structure is highly non- 
monotonic, making the former an ill-suited parameter for describing progen- 
itor structure in core collapse simulations. 

nuclear equation of state. The empirical parameter introduced 
in[0'Connor & Ott (2011) is the compactness of the progen- 
itor, measured at the time of core bounce. It is an inverse 
measure of the radial extent of a given mass coordinate at the 
time of bounce. 



M/Mq 



^(Mbary =M)/1000km 



(12) 



where /?(Mbary = M) is the radial coordinate that encloses a 
baryonic ma ss of M at the time of core bounce. In |Q' Connor 
& Qtt| p011), we chose M = 2.5M0, since this is the rele 



vant mass scale for black hole formation, i.e., a typical max- 
imum baryonic mass at which a range of EOS can no longer 
support a neutron star against gravity. In this study, we pri- 
marily use ^1.75. The motivation for this is that during the 
postbounce preexplosion phase, the relevant mass scale, espe- 
cially for models with relatively small compactness, is much 
less than 2.5 M©. In this study, we choose 1.75 because 
this is close to the average baryonic mass inside the shock at 
200-300 ms after bounce for all models: in the two extreme 
models that span the space in compactness parameter (model 
S12WH07, [^1 75 = 0.24 and ^25 = 0.022], on the lower end; 
model S40WH07 [6.75 = 1-33 and ^2.5 = 0.59] on the upper 
end), the baryonic mass accreted through the shock at 250 ms 
after bounce is 1.45 M and 2.05 M, respective ly. W e further 
justify our motivation of using ^1.75 over ^2.5 in |5.1[ In Fig.[T] 
we plot both ^1.75 and ^2.5 versus ZAMS mass for all 32 con- 
sidered models. ^1.75 is provided in Table [T] for all models. 

For Fig.[T] one notes that while ^1.75 and ^2.5 differ quanti- 
tatively, there is no significant qualitative difference between 
them. The overall trends transcending individual models re- 
main, including the two regions of high compactness near 
22-25 M0 and 35-45 M©. ^1.75 simply provides a more fine- 
grained parameterization at the lower mass scale relevant in 
the first few hundred milliseconds after bounce. Note, how- 
ever, that there are a few models that have similar ^2.5, but 
rather different density structure at small enclosed masses 
and radii and, hence, a different ^1.75. Models ^14WH07 and 
^16WH07 are examples. 

In this study we perform core collapse simulations with 
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each p rogenitor and two EOS. We use the EOS of |Lattimer & 
[Swest y ( 1991 ) with a nuclear incompressibiUty of 220 Me V. 
The LS220 EOS is based on a compressible liquid-drop model 
of the nucleus. Of the publicly available nuclear EOS, the 
LS220 EOS best matches the constraints from n uclear the- 
ory a nd astrophysical observati ons (see Fig. 1 of 'Ott et al.| 
2QTT| and Demorest et al. 2 0101 [Hebeler et al. 2010; Steiner 



et al.|2010i Ozel et al.„2010|). We also e mploy the relativis- 
tic mean field EOS of |Shen et aL] ( |2011| ) that is based on the 



TMl parameter set. It is very different from the LS220 EOS. 
The maximum neutrino-less /3 -equilibrium cold neutron star 
gravitational masses are 2.04 and 2.24 for the LS220 
and HShen EOS, respectively. The radius of a neutrino-less 
/^-equilibrium cold neutron star with a gravitational mass of 
1.4 M0 using the LS220 EOS is 12.7km. For the HShen EOS, 
the corresponding radius is 14.6 km. For details on our par- 
ticular implementation and the treatment of the low-density 
EOS, we refer the reader to O'Connor & Ott (2010, 2011). 
The EOS tables, reader and interpolation routines are avail- 
able from http : / / www .stellar coll apse. org[ 

4. COMPARISON OF NUGRID TO BOLTZMANN TRANSPORT 

Since our implementation of neutrino transport is new and 
approximate, a comparison with published results of full 
Boltzmann neutrino transport is warranted. This will allow 
us to test the ability of our code to reproduce the neutrino lu- 
minosities and spectral properties in the preexplosion phase. 



We compare nuGRlD with the results of ILiebendorf er et a T] 
([2005), a comparison study between two Boltzmann neu- 
trino transport codee^ Th e two codes, A g//g-BOLTZTRAN 
dLiebendorfer et al.|[20Q4| ), and VERTEX dRampp & Janka 
2002'), approach the neutrino transport proW 



|2002r approach the neutrino transport problem in very dif 
Terent ways. Their results compare very well in the New 
Ionian limit, but show significant quantitative differences in 
the general relativistic case. Subsequent modifications to the 
approximate general relativistic potential used in VERTEX 



( |Marek et al.||2005]) have since removed many of the quan- 
titative differences I petween the codes. Th e general relativis- 



tic test case of 'Lieben dorfer et "aL] ( |2005 ) was the collapse 
and early postbounce ev olution of a 15-Mq (a t ZAMS) solar- 
metallicity progenitor of |Woosley & Weaver] ( | 1995 ), referred 
to as model >s'15 WW95 in the following. They employed 
the LSI 80 EOS ( [Lattimer & Swesty] [T^T] ) and a baseline 
set of neutrino-matter interactions, including coupling of en- 
ergy groups via inelastic scattering processes (Bruenn 1985] 
Liebendorfer et al.|2005). We repeat their test here, using the 



same initial conditions and EOS, with our current approxi- 
mations and compare the neutrino observables. We empha- 
size again that the current version of nuGRlD lacks inelastic 
neutrino-electron scattering and velocity depen dent transport 
terms. Both are included in the simulations of "Liebendorfer] 
]et al.j (|2005|). Our transport scheme evolves only the zeroth 
and first moment of the neutrino distribution function, using 
an analy tic closure to truncate the se ries of moment equations, 
whereas ]Liebendorfer et al.) (|2005j solve the full Boltzmann 
equation for neutrino transport. 

In Fig. ]2] we show the luminosities (top panels) and root 
mean squared energies (bottom panels) of three neutrino 

^ The numerical d ata fro m this study are available online at 

\unhbox\voidb @x\hbox{htt p : //iopscience . iop . org/ | 
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Fig. 2. — Neutrino luminosities and root mean squared energies plotted as a 
function of postbounce time for the 515WW95 progenitor. These luminosities 
(top panels) and energies (bottom panels) correspond to the comparison study 
of Liebendorfer et al. ( 2005 1. The left panels contain results for Ue, the right 
panels show h>e (thick lines) and Ux (thin lines) results. The inset plot in the 
upper left panel shows the Ue luminosity around core bounce. Shown in solid 
black lines are luminosities and root mean squared energies obtained with 
nuGRlD. The blue dashed lines and red dashed-dotted lines are the results 
from Liebendorfer et al. (2005 1 using the VERTEX code ( Rampp & Janka] 
12002) and Ag//e-BOLTZTRAN code (Liebendorfer et al. 2004}, respectively. 
A detailed discussion of the differences is provided in the text. 

species: Ve (left panels), and Ve and Vx (right panels). Both 
the luminosity and r oot mean squared energies are defined in 
]Liebendorfer et al.] 12005), Section 4. The black solid lines 
are the results obtained with nuGRlD, the red dashed-dotted 
lines are the predictions of Ag//^-BOLTZTRAN, and the blue 
dashed lines are the VERTEX results. Overall, the agreement 
is good, however there are are several systematic differences: 

(/) The magnitudes of the Ve and luminosities in the early 
postbounce phase predicted by nuGRlD agree well with the 
Ag//^-BOLTZTRAN results but they are systematically lower 
than the VERTEX results. This discrepancy, which also ex- 
ists between VERTEX and the A^//g-BOL TZTRAN, has been 
further investigated inl Marek et al.] (|2005|) and has since been 
resolved. The updated VERTEX code employs an improved 
general relativistic potential and gives comparable amplitudes 
to Ag//^-BOLTZTRAN, and hence nuGRlD. 

(//) The time when the silicon-oxygen interface (located at 
a baryonic mass coordinate of 1.43 in model ^15WW95) 
accretes through the shock, which is marked by a sharp drop 
in the Ve and Ve luminosities, is earlier ( at ^ 140 ms) in our 
simulations than in the simulations of ILiebendorfer et all 
(j2005) {r^ 180ms in the VERTEX simulation^. We have car- 
ried out a number of tests to try to understand this difference. 
We have varied the zero point of the internal energy, which 
affects the relativistic enthalpy entering the momentum and 
energy equations and alters the collapse time. We have at- 
tempted different mappings of the presupernova stellar struc- 

^ Such a sharp d rop is not seen in the A^//^-BOLTZTRAN results. 
[Liebendorf er et al] (|2005j attribute the lack of a sharp drop in the 
Agile- B L T Z T RAN results to the use of an adaptive grid, which introduces 
artificial diffusion and smears out sharp density features. Nevertheless, the 
slow decline of the Ue and Ue luminosities begins around the same time as in 
the VERTEX simulations. 
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ture to our Eulerian grid, and we have tested variations in pre- 
collapse electron capture by par ameterizing as a function 
of density in the collapse phase ( Liebendorf er|20Q5| . None of 
these tests led to a change of the postbounce time to silicon- 
oxygen interface accretion by more than ~ 10 ms. 

Without currently having the modeling technology to test it, 
we suspect that the apparent difference may be due to the only 
other obvious difference between nuGRlD and VERTEX: our 
simplistic treatment of the EOS at low densities. Below a 
density of 6 x 10^ gcm~^, VERTEX replaces the nuclear sta- 
tistical equilibrium (NSE) EOS with an EOS that specifically 
depends on the composition of the matter. Ag//^-BOLT Z TRAN 
assumes all matter below this density is silicon and bums it to 
NSE when the temperature reaches .44 MeV ( Li ebendorf er| 
|et al.|2 005 ). In GRID, as discussed in O'Connor &O tt (2010), 
we assume NSE compositions from the nuclear EOS. At den- 
sities below the validity regime of the nuclear EOS, we take 
the compositions at the lowest density point from the nuclear 
EOS and use the Timmes EOS (Timmes & Arnett 1999 ). For 
reference, the initial density of the silicon-oxygen interface is 
0.4-1. Ox 10^ g cm~^, where the range represents the extent in 
density space. 

(///) The and i>e root mean squared energies, predicted 
by nuGRlD agree very well with the Boltzmann transport re- 
sults during the postbounce phase. The difference seen in the 
i^x root m ean squared energ y is sim ilar to that observed by 
[Thomp son (2002 ) and Lentz et al.|f2012a) when investigat- 
ing the effects of inelastic neutrino-electron scattering. In the 
postbounce evolution, this interaction is expected to predom- 
inately affect the neutrino. We currently ignore this pro- 
cess in nuGRlD and note that the Ux luminosity predicted by 
nuGRlD still agrees well with the full Boltzmann results. 

(iv) Another difference between the evolution in nuGRlD 
and the full Boltzmann transport results arises in the collapse 
phase. The lack of velocity terms and inelastic Ue - e~ scat- 
tering significantly effects the composition of the inner core. 
In simulations with inelastic Ue ~ ^ scattering, neutrinos from 
electron capture on free protons down-scatter off of electrons 
to lower energies. Since the optical depth is lower, these neu- 
trinos can then escape, deleptonizing the core. In our simu- 
lations, these high-energy neutrinos cannot down- scatter and 
therefore cannot escape. Deleptonization is suppressed until 
later phases, when central density and temperature are higher. 
The lack of velocity-dependent terms delays full trapping that 
would normally begin to occur at p > 1 x lO^^gcm"^ until 
nuclear densities, allowing for further deleptonization. At 
bounce the central value of are ~ 0.22 compared to ~0.29 
in VERTEX and Ag//^-BOLTZTRAN. The root mean squared 
Ue energies predicted by nuGRlD are higher during the pre- 
bounce phase, because the neutrinos do not experience the 
down-scattering via inelastic neutrino-electron scattering. 

We find that differences in the inner core structure and com- 
position at bounce do not manifest themselves in the neu- 
trino signal after the collapse phase. Therefore, for this study, 
we find the current version of nuGRlD to be acceptable, 
since our primary focus is the neutrino signal of the preex- 
plosion accretion phase. We do note, however, that the lack 
of energy-coupling terms in our transport can cause quali- 
tative differences near black hole formation. When energy- 
coupling terms, such as gravitational redshift, are included, 
then the and luminosities dr op off at times very close to 
black hole formation, as seen in lFischeF et al., ( ,2009| ). This 



TABLE 1 
Key Neutrino Quantities 



Model 


6.75 


^400ms^^400ms^^400ms 


yyibd 
200ms 


yyibd 
400ms 










LS220/HShen 


LS220/HShen 






[B] 




ri n3i 


ri n^i 


sl2 


0.235 


19.24/14.19/ 


7.73 


1.02/0.92 


2.13/1.78 


sl3 


0.383 


22.18/ 16.58/ 


8.76 


1.25/1.09 


2.53/2.07 


sl4 


0.537 


25.19/ 19.35 / 


9.24 


1.36/1.20 


3.06/2.49 


sl5 


0.580 


25.51 / 19.59/ 


9.17 


1.30/1.16 


3.13/2.59 


sl6 


0.338 


18.91 / 13.72/ 


8.20 


1.11/0.95 


2.00/1.68 


sll 


0.383 


19.93 / 14.54/ 


8.57 


1.20/1.02 


2.13/1.78 


^18 


0.738 


28.66 / 22.26 / 


10.30 


1.55/1.36 


3.62/2.92 


^19 


0.544 


23.36/ 17.49/ 


9.37 


1.43/1.23 


2.67/2.18 


^20 


0.944 


29.39/22.70/ 


11.08 


1.79/1.55 


3.64/2.91 


s21 


0.325 


18.48/13.41 / 


8.09 


1.04/0.91 


1.95/1.65 


s22 


0.972 


30.06/23.29/ 


11.26 


1.82/1.58 


3.76 / 3.00 


s23 


1.256 


42.00/33.88/ 


14.68 


2.29 / 2.00 


5.99/4.67 


s24 


1.167 


39.14/31.35 / 


13.66 


2.16/1.89 


5.44/4.24 


s25 


1.040 


32.30/25.28/ 


11.87 


1.93/1.68 


4.15/3.28 


s26 


0.727 


24.97/ 18.84/ 


10.10 


1.59/1.38 


2.88/2.33 


s21 


0.783 


25.67 / 19.45 / 


10.21 


1.61/1.40 


3.01/2.42 


^28 


0.640 


25.86/ 19.83/ 


9.54 


1.36/1.21 


3.16/2.62 


s29 


0.556 


23.24/ 17.38/ 


9.43 


1.45/1.22 


2.64/2.17 


^30 


0.760 


26.61 / 20.30/ 


10.19 


1.60/1.40 


3.20/2.58 


^31 


0.687 


25.23 / 19.09 / 


9.95 


1.56/1.35 


2.96 / 2.40 


^32 


0.883 


28.10/21.58/ 


10.70 


1.71/1.49 


3.43/2.75 


^33 


0.965 


30.36/23.57/ 


11.28 


1.82/1.58 


3.82/3.04 


^35 


1.129 


35.93/28.50/ 


12.77 


2.06/1.79 


4.83/3.79 


^40 


1.328 


48.32 / 39.45 / 


17.51 


2.49/2.21 


7.23/5.75 


^45 


1.300 


46.42/37.81 / 


16.56 


2.42/2.14 


6.86/5.44 


^50 


0.701 


25.60/ 19.42/ 


10.00 


1.57/1.32 


3.02/2.45 


^55 


0.577 


23.56/ 17.65/ 


9.52 


1.48/1.25 


2.69 / 2.20 


^60 


0.461 


21.52/ 15.88/ 


9.07 


1.33/1.12 


2.36/1.95 


^70 


0.755 


25.83/ 19.58/ 


10.17 


1.60/1.39 


3.04/2.46 


580 


0.591 


23.75/ 17.80/ 


9.63 


1.48/1.26 


2.71/2.22 


5100 


0.792 


29.90/23.33/ 


10.59 


1.61/1.42 


3.82/3.07 


5120 


0.474 


22.31 / 16.64/ 


9.02 


1.33/1.15 


2.52 / 2.07 



Note. — For each model in the Woosley & Heger (2007) model set 
we show 6.75, the cumulative emitted neutrino energy in z/g, and a 
single Ux at 400 ms after bounce. The numbers correspond to the models 
run with the LS220 EOS. We also present, for each model, the estimated 
number of IBD events in a Super-Kamiokande-like detector at 200 and 
400 ms after bounce for a supernova at a fiducial galactic distance of 
lOkpcfor both EOS. 

effect can be captured with GRlD's l eakage scheme in which 
it is trivial to include redshift terms ([ O'Connor & Qtt|201Ql ). 
Capturing redshift in an energy-dependent transport scheme 
requires energy-group coupling. We have also compared our 
Ml scheme to the results of Fischer et al. ( 2009 ) for the 40 
model from Woosley & Weaver] ( |1 995 j ) using the LSI 80 EOS. 
In this model a black hole forms within 500 ms of bounce. We 
find differences in the Ve and Ve luminosities of ~ 10-20% in 
the last ~50ms, but good agreement in the early postbounce 
phase. 

5. RESULTS 

5.1. Trends in the Neutrino Observables 

We perform core collapse and early postbounce evolutions 
of all 32 models introduced in fusing both the LS220 and 
the HShen EOS. In Fig. [3] we present the three neutrino lu- 
minosities and average energies for each model and EOS as a 
function of postbounce time. We do not expect a clear trend in 
the neutrino observables with ZAMS mass. However, we do 
expect trends based on the presupernova structure of the star, 
which is well encapsulated by the compactness parameter in- 
troduced in ^ In the top set of panels we show simulations 
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run using the LS220 EOS. In the bottom panels, simulations 
performed with the HShen EOS are shown. To highlight that 
there is indeed a trend with presupernova structure, we color- 
code individual models according to their compactness pa- 
rameter. The mapping between line color and ^1.75 is provided 
on the right. To more directly highlight the EOS dependence, 
we include the luminosity and average energies of two models 
run with the LS220 EOS in the HShen EOS panels with thick 
dashed lines. These models, ^12WH07 and ^40WH07, have 
the lowest and highest compactness parameter in our model 
set, respectively. 

We find that there is little variation in the peak luminosity of 
the Ve neutronization burst signal. For all 32 models simulated 
using the LS220 (HShen) EOS, the peak amplitude varies by 
less than 3% (5%) from the average. This reflects the univer- 



sal nature of the collapse of the inner core ( [Liebendorfer et al. 
p002). After the neutronization burst, the postbounce lumi- 
nosities of all species increase systematically with increasing 
compactness parameter. Models with higher ^1.7 5 have higher 
temperatures throughout the protoneutron star ( [O'Connor &| 
[Ott 2Qllj . This increases the diffusive neutrino luminosity 
and" is best seen in the Vj^ luminosities. The postbounce ac- 
cretion rate also increases with the compactness parameter. 
Higher accretion rates, and the deeper gravitational poten- 
tial due to the higher protoneutron star mass, increase the ac- 
cretion luminosity, which is most directly reflected in the 
and Ve signals. After roughly 100 ms, the average energy of 
the emitted neutrinos also shows an increasing trend with the 
compactness parameter. The matter temperature at the neutri- 
nosphere is higher in models with larger compactness, there- 
fore a higher average neutrino energy is observed at infinity. 

The neutrino luminosities and average energies from sim- 
ulations using the HShen EOS are systematically lower than 
the luminosities and average energies from simulations of the 
same model run with the LS220 EOS. This is clearly seen in 
the bottom set of panels in Fig. [3] models ^12WH07-LS220 
and ^40WH07-LS220 have luminosities and average energies 
that are comparable to or larger than in the corresponding 
HShen models. For a fixed accretion rate (or fixed progeni- 
tor model), the location of the neutrinosphere of each species 
influences the emitted luminosities and spectra. In models 
evolved with the stiff HShen EOS, the neutrinospheres are 
located systematically at larger radii and lower matter tem- 
peratures than in models run with the softer LS220 EOS. For 
example, in the ^12WH07 simulations, the Rosseland-mean 
Ve neutrinospheres at 200 ms after bounce have radii and tem- 
peratures of ~ 35.3 km and ~4.86MeV; and ~ 39.5 km and 
- 4.46 MeV for the LS220 and the HShen EOS, respectively. 
The larger neutrinosphere radii are responsible for the lower 
accretion luminosity since the latter is set essentially by the 
product of the mass accretion rate and the gravitational po- 
tential at the protoneutron star surface (Liebendorfer et al.| 
[2002 ). The latter is located at larger radii in simulations us- 
ing the HShen EOS. The differences in the neutrinosphere 
radii and temperatures between the LS220 and HShen EOS 
also give an explanation for the systematically lower aver- 
age neutrino energies seen in the HShen simulations. Mat- 
ter at larger radii has been compressed less and therefore is 
cooler. This leads to average neutrino energies that can be 
up to 5 MeV lower for the HShen EOS than for the LS220 
EOS for the same progenitor model (see Fig. |3]). The dif- 
ference in the neutrino luminosities and average energies be- 
tween the two EOS is largest for models with ^1.75 > 1.2, in 



which the high accretion rates lead to the accumulation of 
^2Mq of material inside the shock within ^ 200-300 ms of 
bounce. In the case of the LS220 EOS, this leads to very high 
temperatures throughout the protoneutron star as it becomes 
more and more compact and closer to gravitational collapse 
to a black hole. In our simulations, the most compact model 
^40WH07 (^1.75 = 1.33, ^2.5 = 0.59) forms a black hole 503 ms 
after bounce. The slightly less compact model ^45WH07 
(6.75 = 1-30, 6.5 = 0.55) forms a black hole 563 ms after 
bounce. The high temperatures present in the LS220 simu- 
lations at these times will not be obtained until postbounce 
times > 1 s in models using the HShen EOSj^ 

Most neutrino detectors are most sensitive to the electron 
antineutrino luminosity through the dominant IBD interac- 
tion, Ve+p ^ n-\-e^. In the left panel of Fig. |4j we consider 
the cumulative emitted Ve energy for each model using the 
LS220 EOS. We color code the models based on their com- 
pactness parameter and include two reference models that use 
the HShen EOS, model ^12WH07 and model ^40WH07. The 
graphs shown in this panel are the integral of the graphs shown 
in the top center panel of Fig. [3] It is obvious that the cumu- 
lative amount of Ve energy emitted during the preexplosion 
phase strongly correlates with the compactness of the progen- 
itor model. For example, the amount of emitted Ve energy 
from model ^40WH07 (6.75 = 1.33) is always between two 
and three times of that of model ^12WH07 (6.75 = 0.24). We 
make this point more quantitative in the center and right pan- 
els of Fig. |4] In the center (right) panel we plot the cumu- 
lative emitted Ve energy at 100, 200, 300, and 400 ms after 
bounce for both EOS as a function of 6.75 (6.5). For refer- 
ence, we present a subset of these numbers in Table [T] We 
see a very clear correlation that depends only weakly on the 
chosen EOS. Note, however, that for models with small com- 
pactness parameter (6.75 < 0.8), the correlation between the 
total emitted Ve energy and the compactness parameter is not 
as strong after 400 ms of postbounce evolution. Comparing 
the center and right panels justifies our choice of 6.75 over 
6.5 as explained in ^ 

The onset of an explosion will break the correlation ob- 
served in Fig. [4] Once it is launched, the accretion lumi- 
nosity effectively turns off and only the diffusion luminosity 
remains. One also expects this diffusion luminosity to show 
a correlation with the compactness of the progenitor, since 
the remnant protoneutron star's thermodynamic conditions, 
such as the central entropy and its mass, are essentially set 
by the presupernova structure. However, it is currently un- 
clear whether one should obtain a correlation between explo- 
sion time and the compactness parameter. Clarification will 
require a more complete understanding of the core-collapse 
supernova explosion mechanism and may require extensive 
parameter studies with fully self-consistent three-dimensional 
simulations. 

5.2. Detectability 

We use the publicly available software SNOwGLoBES Qto 
predict the neutrino signal observed in Earth-based neutrino 

^ While we d o not follow these mode ls to black hole formation in our 
current study, in [O'Connor & Ott| (|2011| we found that the black hole for- 
mation times of the ^4U WH07 and k-5WH07 models are ~1.3 s and ~1.4s, 

respectively, when using the HShen EOS. 

available at http : / /www . phy . duke . edu/~schol/| 
snowglobes 
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Fig. 3. — Neutrino luminosities (top panels) and average energies (bottom panels) plotted as a function of postbounce time for all 32 models of |Woosleyl 
|& Heger ( 2007). The top set of panels shows results obtained with the LS220 EOS. The bottom panel shows the same for the HShen EOS, but includes, tor 
reference, two LS220 models: 5'12WH07 and 5'40WH07. The left, center, and right panels show results for Ue, i^e, and Ux, respectively. The curves are color- 
and line- weight-coded with increasing compactness (^1.75), the mapping from color to compactness parameter is shown on the right. There is a clear trend in 
all luminosities and average energies with compactness parameter. The progenitor with the highest compactness, 540WH07, forms a black hole at 502 ms after 
bounce. None of these models explode, but the onset of an explosion in any of these models would lead to a sudden deep drop (strongest for Ue and i?e) in the 
luminosities and average energies (F ischer et al. 20 10). The smaller drop observed for most models models here is due to the sudden decrease of the accretion 
rate when the silicon-oxygen interface reaches the stalled shock. 



detectors. SNOwGLoBES (Scholberg'2012), which in turn re- 
lies on GLOBES (Ruber et al. 2005, 2007), is a set of rou- 
tines that compute the interaction rates of supernova neutri- 
nos in user- specified detector configurations. Variables in- 
clude detector material (e.g., water, scintillator, argon, and 
lead), detector volume, detector response functions, and a 
host of relevant neutrino interactions. For this investigation, 
we consider only IBD events in a water Cherenkov detec- 
tor. We choose a detector mass of 32 kT, the mass of wa- 
ter in Super-Kamio kande that is sens itive to core-collapse su- 
pernova neutrinos (Scholberg 2012| ). For reference, we use 
the wcl00kt30prct smearmg rates and efficiencies pro- 
vided with SNOwGLoBES. We construct SNOwGLoBES ini- 
tial fluence data from our simulations binned in 5 ms inter- 
vals. We provide these energy-dependent fluences at http : 
//www . stellarcollapse . org/Mlprog for all mod- 
els, neutrino species, and both EOS in 5 ms intervals up to 
450 ms after bounce. We assume a fiducial galactic supernova 
distance of lOkpc. In Fig.|5] we show the cumulative number 



of events for each model run with the LS220 EOS and, for 
reference, two models run with the HShen EOS. The lines are 
color-coded according to compactness parameter. Note that 
the vertical scale in this figure is in thousands of events. We 
note that our answers agree with the total number of expected 
events in Super-Kamiokande from a galactic core-c ollapse su- 
pernova at 10 kpc, which is estimated to be ~ 7000 ( [Scholberg 
2012). To arrive at this number from our results, consider 
the lowest ZAMS mass progenitor in our model set, model 
^12WH07. After 450 ms of evolution, 15B of energy has 
been radiated (Fig.|4]), which corresponds to 2000 IBD events 
(Fig. 0). For 50 B of released energy (- 1/6 of 300 B, the 
fiducial energy released in neutrinos over the entire cooling 
phase), one would then expect ~ 7000 events. However, as is 
clear from Fig. |5j the number of events from the next galac- 
tic supernova may be higher than this fiducial number. More 
importantly, the rate of events in the preexplosion phase will 
give us detailed information on the progenitor core structure. 
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Fig. 4. — Cumulative e mitted electron antineutrin o energy. We show the cumulative totals via several methods. In the left panel we show time series data for 
each of the 32 models of'Woosley & Hegerl (120071 run with the LS220 EOS. The color coding corresponds to the value of ^1.75 and is given in Fig. [3] This 
cumulative energy is the time integral of the top-center panel of Fig. [3] The dashed black lines correspond to models 512WH07 and ^40WH07 run with the 
HShen EOS. In the center and right panels we present the emitted energy at select postbounce times, for each model and EOS, plotted versus the compactness 
^1.75 (center panel) and ^2.5 (right panel) of the model. 

on the detector's response function and efficiency. 

An additional independent path to experimentally probing 
the inner structure of the progenitor is via the total neutrino 
energy emitted in all species over the first 10s of seconds 
after the initial collapse. This method requires a measure- 
ment of the total fluence of neutrinos of all species, not just 
electron antineutrinos. An example of a neutrino interaction 
capable of relaying such information is the mono-energetic 
de-excitation of a neutral-current neutrino-excitation of ^^C 
( [Scholberg 2012 ). Such a measurement would require good 
energy resolution and a significant source of carbon, for ex- 
ample, a liquid scintillator neutrino detector. We note that 
even with a liquid scintillator detector, the dominant neutrino 
interaction is still IBD ( |Scholberg|2012| ). 

If such a measurement was made, and there is not a sig- 
nificant amount of rotation (see the discussion on rotation in 
|5.3| ), one can immediately infer the gravitational binding en- 
ergy of the remnant, since neutrinos carry away the vast ma- 
jority (~ 99^J) of the gravitational binding energy. For typ- 
ical nuclear EOS like the ones considered here, this results 
in a one-to-one mapping of the released gravitational binding 
energy to the baryonic mass of the remnant, and, hence, the 
gravitational mass of the remnant. This is most easily seen by 
fitting the gravitational binding energy of a cold (T=0. 1 MeV), 
neutrino-less /^-equilibrium, non-rotating neutron star to its 
baryonic mass. From cold neutron star TOV solutions using 
the LS220 EOS one can obtain an empirical fit to better than 
3% above a baryonic mass of 1.15 M©, 




100 



200 300 
t-tbounce [n^s] 



400 



Fig. 5. — Cumulative IBD events in a Super- Kamiokande-like water 
Cherenkov detector at a fiducial galactic distance of lOkpc versus postbounce 
time. We use the SNOwGLoBES package to determine the integrated IBD 
event rate in a 32 kT water Cherenkov detector at lOkpc. The color coding 
corresponds to the value of ^1 75 and is provided in Fig.l3] The dashed lines 
are results for models ^12WH07 and ^40WH07 run with the HShen EOS. 
In the inset we show the cumulative IBD hits as a function of ^1.75 for each 
model and EOS at four postbounce times: 100, 200, 300, and 400 ms. 

In order to more directly quantify the differences between 
variations in progenitor compactness and variations in the 
nuclear EOS, we plot in the inset of Fig. [5] the number 
of expected IBD events in a Super-Kamiokande-like water 
Cherenkov detector at various postbounce times versus ^1.75. 
There is a well defined trend: The number of IBD events de- 
tected in the first 100, 200, 300, and 400 ms increases with 
the compactness parameter of the models. For reference, we 
include the expected number of events at 200 and 400 ms for 
both EOS in Table [T] We find that the EOS dependence of 
the expected number of events is similar to the EOS depen- 
dence of the total emitted i>e energy: the HShen EOS leads 
to a lower number of events (compare the inset of Fig. |5] to 
the center panel of Fig.|4]). The dependence on EOS is some- 
what stronger here, since the the lower average energy pre- 
dicted from stiffer EOS translates into a reduced cross section 
in Earth-based detectors. In addition to the total number of 
events, a water Cherenkov detector measures individual ener- 
gies, and thus, allows for the reconstruction of the cumulative 
emitted Ue energy over time. This reconstruction will depend 



1.12 xlO^^(Mbary/M0)2 ergs. 



(13) 



^binding ~ i.iZ X iU V^Kibary 

A similar fit for the HShen EOS gives, 

^bmdmg - 9.78 X lO^^^Mbary /Mo)^ ergs , (14) 

and is accurate to 5% above baryonic masses of 1.15 Mq. Be- 
low Mbary = 1.15M0, the empirical quadratic fit is not as ac- 
curate. However, all models considered here reach a baryonic 
protoneutron star mass of 1.15 within ~ 10 ms of bounce. 
Hence, we believe that the above fits are acceptable for the 
iron-core core collapse events considered here. 

^ The remaining ~1% of the energy is predominantly shared among the 
kinetic energy of the explosion, the original binding energy of the unbound 
stellar mantle, and the binding energy of the iron core at the onset of collapse. 
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We now make the assumption that an explosion launched 
at a particular postbounce time will result in a neutron star 
remnant with a baryonic mass equal to the baryonic mass that 
has accreted through the shock up until the time of the ex- 
plosion. This neglects any late-time fallback of material onto 
the protoneutron star which would lead to additional neutrino 
emission. Fryer ( 2009 ) and Ugliano et al.| ( |2012] ) predict fall- 
back masses <5-10% of the initial protoneutron star remnant 
mass. We also neglect any asymmetric mass accretion which 
may occur in the early explosion phase. In Fig.|6] we convert 
the baryonic mass enclosed by the shock to the total emitted 
neutrino energy using Equations [13] and [14] We plot this for 
all progenitor models (run with the LS220 EOS) as a function 
of the hypothetical time of explosion. As a concrete exam- 
ple, consider the situation where 300 B (shown as the dashed 
line in Fig. [6]) of total neutrino energy was observational in- 
ferred. This could correspond (/) to a progenitor with a high 
compactness parameter that exploded at an early time, e.g., 
model ^40WH07 at 70 ms or (//) to a low progenitor with 
low compactness parameter that exploded at late times, e.g., 
model ^12WH07 at 400 ms. If we have an estimate of the ex- 
plosion time, e.g., via characteristic features in the neutrino 
observables, then we can use the combined measurement to 
probe the progenitor core structure. This is further quantified 
in Fig.^ where we choose three total emitted neutrino ener- 
gies, 250, 300, and 350 B, and determine the time at which the 
explosion must have been launched for a given compactness 
and the respective total emitted energy. We plot this explo- 
sion time versus ^1.75 for all models and both EOS. In gen- 
eral, for a fixed total emitted energy, as the compactness pa- 
rameter of the progenitor increases, the explosion time must 
decrease. For progenitors with high ^1.75 (> 0.8), there is a 
clear mapping between the explosion time and the compact- 
ness parameter, given a specific total emitted neutrino energy 
and an EOS. As was the case for the total IBD rates, there is 
some scatter at low ^1.75 (< 0.8), but there is still an overall 
trend. 

Fig. [T] shows that there is a very strong EOS dependence. 
While the baryonic mass inside the shock as a function of 
postbounce time does not vary strongly with EOS, the grav- 
itational binding energy released does. Based on the above 
empirical fits to the gravitational binding energy, the LS220 
EOS leads to ~ 14% more energy release than the HShen 
EOS for the same baryonic mass. This method of determin- 
ing the compactness parameter (by combining an estimate of 
the explosion time and the total emitted neutrino energy) can 
be used together with the methods described above (using the 
postbounce, preexplosion IBD rates) as a consistency check, 
or to break degeneracies, which we will discuss next in |5.3| 

5.3. Degeneracies in Neutrino Observables 

There are a number of degeneracies and uncertainties that 
may prevent fully conclusive statements regarding the map- 
ping from detected signal to progenitor core structure. These 
include, nuclear EOS, rotation, viewing angle, distance, and 
neutrino oscillations (including collective oscillations). We 
will discuss each one of these consecutively and indepen- 
dently, although all may be relevant in a generic situation. 

Nuclear EOS: We have already briefly explored the depen- 
dence of the neutrino observables on the nuclear EOS by tak- 
ing two very different EOS and comparing the emitted neu- 
trino signal. The total number of events predicted to be de- 
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Fig. 6. — Total emitted neutrino energy estimated from the enclosed bary- 
onic mass as a function of explosion time for the LS220 EOS. A mea- 
surement of the total emitted neutrino energy and an estimate of the ex- 
plosion time constrains progenitor structure. For reference, we provide the 
cold neutron star gravitational mass associated with the released binding 
energy on the right ordinate. This figure is constructed using a fit of the 
gravitational binding energy of a cold neutron star to its baryonic mass, 
E^S^^^ ~ 1.12 X lO^^(Mbary/>W0)^ergs, and the baryonic mass enclosed in 
the shock at any given time. This defines the explosion time to be, in a La- 
grangian sense, the time at which the outermost final neutron star mass ele- 
ment accretes through the shock. The color coding corresponds to ^1.75, the 
color coding is provided in Fig.|3] 
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Fig. 7. — Explosion time for various progenitors and EOS, under the con- 
straint of a fixed amount of released neutrino energy. We show three values 
of the total released energy, 250, 300, and 350 B. Progenitors with high com- 
pactness must explode at earlier times to achieve the same total amount of 
released energy in neutrinos when compared to a progenitor with a lower 
compactness. The strong EOS dependence is due to the different binding en- 
ergies resulting from the two EOS. The LS220 EOS leads to a more compact 
neutron star, releasing ~ 14% energy for a given baryonic mass. 



tected from a given progenitor model varies with EOS (cf., 
Fig. [5|) in such a way that a high compactness model paired 
with the HShen EOS produces the same number of events as a 
model with a slightly lower compactness paired with a softer 
EOS, like the LS220. Hence, there is a clear degeneracy be- 
tween the progenitor compactness and the equation of state. 
However, a water Cherenkov detector will not only detect a 
given number of IBD events, but also their energy distribu- 
tion. From this distribution, and the detector response, one 
can work backwards to reconstruct both the emitted average 
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Fig. 8. — Average i>e energy versus cumulative emitted i>e energy each EOS 
and progenitor. This is shown for three postbounce times, 100 ms (left panel); 
200 ms (center panel); and 300 ms (right panel). Solid (open) points corre- 
spond to the LS220 (HShen) EOS. The degeneracy between EOS and cumu- 
lative Ue energy is broken when the average energy is taken into account. We 
find that is a robust prediction across all progenitors and postbounce preex- 
plosion times. 

Ve energy and the luminosity. We show in Fig. [8] how knowl- 
edge of the emitted energy spectrum can break the degeneracy 
with cumulative emitted De energy and EOS, provided the dis- 
tance to the core collapse event is known (see the discussion 
of distance uncertainties later in this section). For each model 
and EOS, we plot the average Ve energy as a function of cu- 
mulative emitted Ve energy at three select postbounce times: 
100, 200, and 300 ms. The data from each of the 32 progeni- 
tors fall on a unique line that is parameterized by the compact- 
ness. Provided large enough statistics, this allows to break the 
EOS-progenitor degeneracy. 

Rotation and Viewing Angle: The effect of rotation on the 
neutrino signal may be more difficult to disentangle from 
the effects of progenitor structure. To explore the effect of 
rotation, we perform 1.5D simulations of rotating core col- 
lapse. These simulations treat rotatio n in a spherically sym 
metric way, using 'shell ular rotation' ( [Thompson et al.|2005 



[O'Connor & Ott 2010| ). This approximation only captures 



the spherically-averaged centrifugal effect of rotation on the 
matter. As the amount of rotation is increased, the protoneu- 
tron star becomes more and more centrifugally supported. 
Increased rotation leads to lower densities and temperatures 
throughout the protoneutron star. This in turn effects the neu- 
trino signal. In Fig. |9] we plot the De luminosities and av- 
erage Ve energies determined from nuGRlD for 12 rotating 
core collapse simulations. We use model ^15WH07 paired 
with the LS220 EOS. T he initial rotation rate is assigned via 
( [O'Connor &Ott|2011[ ) 



j(r) = j 



16, cx) 



1 + 



X 10 



cm' s- 



(15) 



where Amq is the radius that encloses 1 M^t) , which for the 



0' 



^15WH07 progenitor is 703 km, and 7i6,oo is the specific 
angular momentum at infinity, in units of 10^^ cm^ s~^ 
We vary ji6,oo from to 3 in increments of 0.25. For 
the ^15WH07 progenitor, the initial central angular veloc- 
ity is = 2.03 Ji6,cx) rads~\ giving a range of Qc from to 
6.08 rad s~^ . At bounce, the 7i6,oo = 1 , = 2.03 rad s~^ model 
has a rotation period of ~ 2.4 ms, which decreases nearly lin- 
early with increasing initial 7i6,oo- Increased rotation leads 
to the overall reduction of both the luminosity and av- 
erage energy. For example, with Qc = 0, 2.03, 4.05, and 
6.08 rad s~^ the z/^ luminosity at 100 ms is ~55, ~52, ~40, 
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Fig. 9. — Luminosities (top panel) and average energies (bottom panel) 
from 1.5D rotating core collapse simulations of model ^15WH07 for vari- 
ous initial central angular velocities. We assign angular momentum to the 
^15WH07 progenitor via EquationfTslwith values of ji6,oo ranging from to 

3 in 0.25 increments. This gives inrnal central rotation rates of to ~ 6.08 
rad s~^ in increments of ~ 0.51 rad s~^ . In our simulations, increased rotation 
lowers both the neutrino luminosity and the average energy. 

and ~25Bs~\ respectively. While we do not show the Ue 
and Ux luminosities and average energies, they follow the Ue 
trends. The only exception to this is the Ug neutronization 
burst where rotation in our 1.5D approach does not signif- 
icantly alter the luminosity. Our 1.5D approach to rota- 
tion does not capture the anisotropy of the n eutrino emission, 
rather just the angle- averaged value. [Ott et al. ( 2008) consid- 
ered a model with an initial rotation rate of tt rad s~\ roughly 
corresponding to our ji6,oo = 1.5 model. They found that the 
ratio of polar to equatorial luminosity can be as large as 3- 

4 and that average neutrino energies in polar regions can be 
harder by > 1-2 MeV than on the equator. We find that it is not 
possible with the cumulative neutrino signals alone to break 
the degeneracy between rotation and the progenitor star-EOS 
combinations. The angle- averaged Ue luminosity and aver- 
age energy from a rotating system mimic the luminosity 
and average energy from a non-rotating system with a pro- 
genitor with lower compactness and/or a stiff er EOS. In addi- 
tion, the rotational energy stored in the neutron star remnant 
can lead to an underestimate of the total emitted neutrino en- 
ergy, again mimicking a lower-compactness progenitor and/or 
a stiff er EOS. However, all hope is not lost: stell ar evolution 
and current constraints on pulsar birth spins (e.g., [Heger et al.[ 
2005 [ Ott e t al. 2006 ) suggest that most massive stars will be 
spinning too slowly for rotation to have a strong effect on dy- 
namics and neutrino signal. Heger et al.[(l2005 ), for example, 
predict a presupernova rotation rate of 0.2 rad s~^ for a 15 Mq 
star, which is smaller than all rotation rates considered here. 
If, on the other hand, the progenitor is rapidly spinning, the 
degeneracy in the neutrino signal may be broken by coinci- 
dent observations in gravitational waves th at are able to co n- 
strain the rotation rate of the collapse core ( [Ott et al.|2012| ). 

Distance: If the electromagnetic signal is blocked, by, for 
example, the galactic center, we may not be able to obtain 
a reliable distance to the supernova. The Ue flux at Earth 
follows an inverse square law, while the spectral distribution 
will not change over the distances of relevance. This means 
that degeneracy with distance is also hard to break with the 
IBD neutrino signal alone for the following reason: A high- 
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compactness progenitor at a large distance can produce the 
same energy-integrated flux at Earth as a low-compactness 
progenitor. Less neutrinos are emitted in the latter but the 
flux has not been diluted as much due to the closer distance. 
If the nuclear EOS is not known, the observed neutrino flux 
and energy spectra could be associated with a range of nuclear 
EOS-distance combinations. A progenitor at a large distance 
could produce the same neutrino flux and energy distribution 
as a lower-compactness progenitor with a softer EOS. Here, 
the softer EOS, which will give rise to an increased average 
energy, compensates for the decrease in the average energy 
produced by the smaller compactness. 

Future constraints on the nuclear EOS could break this de- 
generacy since one could more reliably associate an observed 
energy distribution with a particular progenitor compactness 
(Fig. [5]). This may even provide a distance e stimate, as would 
a dete ction of the neutronization burst signal ( [KachelrieB et al. 
2QQ5| ). It does not show a strong dependence on the progeni- 
tor model, but is difficult to detect since the cross sections 
in water Cherenkov detectors are much lower than the IBD 
cross section, resulting in few events. Some detector mate- 
rials have significantly larger cross sections and would be 
better suited to detect the neutronization burst, such as liquid 
argon ( Scholberg 2012) and lead (Duba et al. 2008 ). 

Neutrino Oscillations: First, considering only matter- 
induced neutrino oscillations, in the inverted mass hierarchy, 
the electron antineutrino signal at Earth is a composite spec- 
trum of cos^(^i2) ^70% of the original Ue neutrinos and 
sin^(6>i2) ~30% of the original spectrum, where On is 
the mixing angle between the mass eigenstates 1 and 2. In 
the normal neutrino mass hierarchy, the entire Ve signal is 
replaced with the original = Vx signal (Dighe & Smimov 
[2000 ). In either case, we still expect the total number of IBD 
events to increase with the compactness parameter of the pro- 
genitor since both the De and the Ux luminosity increase. If the 
hierarchy of the neutrino mass eigenstates is not determined 
by neutrino experiments before the next nearby core-collapse 
supernova, the early postbounce, preexplo sion neut rino lumi- 
nosities may pro vide an answer ( KachelrieB et al. 12005, Ser-^ 
pico et al^|2012| ). This relies on the systematically different 



rise times of the Vx and Ve signals. 

Much more troublesome are the collective neutrino oscilla- 
tions that arise from coherent neutrino-neutrino forward scat- 
tering. Collective oscillations are very sensitive to the energy 
spectra (both the distribution and magnitude) of all neutrino 
flavors and the background matter density. Since the gov- 
erning equations are highly non-linear, there are currently no 
simple analytic expressions predicting the neutrino signal at 
Earth based on the output of core-collapse simulations. Re- 
cent studies suggest that during the early postbounce, pre- 
explosion phas e, collective n eutrino oscill ations ma y be sup- 
pressed (Chakraborty e t al. 201 1b a; Sarikas etal. 12012), how- 
ever the co mmunity has not yet reach ed consensus, see, e.g., 
[Cherry et al., ( .2012^ and [DasguptTeTal . , ( ,20 1 2^ . 

6. DISCUSSION 

The next nearby core-collapse supernova will be ex- 
tremely well observed in neutrinos. Super- Kamiokande alone 
will observe ~7000 electron antineutrinos from a typical 
core-collapse supernova at a fiducial galactic distance of 
lOkpc. Future detectors of the scale of the proposed Hyper- 



Kamiokande may see in excess of 10^ events. Such high- 
statistics observations will provide rich information on the 
neutrino signal. Comparison with theoretical model pre- 
dictions will allow to falsify or constrain a broad range of 
hypotheses in core-collapse supernova astrophysics and nu- 
clear/neutrino physics. Unexpected signal features may lead 
to the discovery of new physics. 

In this study, our focus has been on the imprint of the 
progenitor star's structure on the neutrino signal in the 
postbounce preexplosion phase of core-collapse supernovae. 
We have carried out a large set of spherically- symmetric 
radiation-hydrodynamics simulations of core collapse and the 
early postbounce phase with the goal of studying trends in the 
neutrino signal with variations in progenitor structure. 

Our re sults show, in agreement with previous work (e. 
Burrows & Mazurek 1983 ; Liebendorfer et al. 2001 2002 



20031 [2 004; Thomps on et al. 2003; KachelrieB et al., 2005 
Buraset al. 2006; Se rpico et ar.|2012l) , that the Ue signal from 
the neutronization burst emerging shortly after bounce has 
very little progenitor dependence, due to the universal nature 
of homologous inner core collapse. 

The neutrino signal in the postbounce preexplosion phase is 
determined primarily by the accretion luminosity of outer iron 
core, silicon shell, and oxygen shell material. The postbounce 
accretion rate depends on the inner structure of the progenitor 
at the presupernova stage. Our results show that preexplo- 
sion neutrino signal has an essentially monotonic dependence 
on progenitor structure described by a single parameter, the 
compactness M/R(M) (where M is a typical baryonic 
mass reaching the center over the timescale of interest). The 
greater a progenitor's the higher are the emitted luminosi- 
ties and average energies of all neutrino species. Scaling in 
the same way is the total emitted energy in neutrinos over the 
entire protoneutron star cooling phase for a given explosion 
time at which accretion is shut off. These trends are robust 
and independent of the nuclear EOS. They are also rather in- 
sensitive to the particular choice of the reference mass M in 
as long as it is in the range of typical neutron star baryonic 
masses (~1.4-2.5M0) and we find ^1.75 to be a good choice. 

The monotonic dependence of the preexplosion neutrino 
emission on progenitor compactness translates directly to the 
neutrino signal observed by detectors, provided collective 
neutrino oscillations do not lead to complicated swaps of fla- 
vor spectra that brake the dependence of the observed sig- 
nal on progenitor structure. Neutrino observations of the next 
nearby core collapse event thus may, in principle, allow quan- 
titative constraints on the inner structure of the progenitor 
star. As an example with real neutrino data, we consider the 
early postbounce neutrino signal o bserved from SN 19 87A 
by the Kamiokande-II experiment (|Hirata et al.||1987|). Of 
the eleven events that were observed, the first four events oc- 
curred within 323 ms of each other. All events obser ved by 
Kamiokand e-II are consistent with being IBD events ( [Hirata 
|et al.|1987] ). We assume that the first event occurs at the on- 
set of the postbounce phase, although the actual bounce time 
is likely somewhat earlier. In Fig.[TO| we plot the cumulative 
number of events observed from SN 1987 A in the first 500 ms 
along with the SNOwGLoBES prediction from our simulations 
with the LS220 EOS using a 2.14kT water Cherenk ov detec- 
tor at 51.47kpc. We use the efficiencies quoted in Burrows 
( |1988| ) and the smearing matrices from the SNOwGLoBES 
detector configuration wclOOktSOprct. To show the full 
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Fig. 10. — Predictions of the cumulative IBD hits in Kamiokande-II for 
a core-collapse event at 51.47 kpc. The left panel assumes no neutrino os- 
cillations, while the right panel assumes a swapping of the Ve and Vx spec- 
tra. The color coded lines denote different progenitors and are generated 
with SNOwGLoBES based on our simulations with the LS220 EOS. More 
details are provided in the text. We overlay the cumulative hits observed in 
Kamiokande-II from SN 1987 A, assuming the first event denotes the time of 
core-bounce. 

range of the possible effects of MSW neutrino oscillations, we 
show the expected number of events assuming no oscillations 
and assuming a complete switch of the and spectra (as 
would be the case in the normal mass hierarchy; the inverted 
hierarchy is a combination of these two signals). While the 
quantitative results obviously depend on neutrino oscillation 
details, the qualitative trend with is unbroken. Compar- 
ing our predictions with the events observed from SN 1987 A 
one notes (but must keep the very small-number statistics in 
mind) that either the explosion must have occurred early in the 
postbounce phase and/or the progenitor must have had a rela- 
tively lo w ^1.75 , which is broadly consistent with the previous 
work of |Bruenn| (|1987 ) and Burrows ( 1988 ). It is also not in- 
consistent with the propos ed ZAMS mass of ~18-2OM0 for 
Sanduleak -69° 202 (e.g., |Woosley etal. 12002] ), the blue su- 
pergiant progenitor star of SN 19 87 A. However, the mapping 
between ZAMS mass and stellar structure (i.e., compactness) 
at the presuper nova stag e appears to be highly n on-monotonic 
(cf. Fig.[TJ Woo sley et a l. 2002 ; Woosley & Heger 2007). This 
makes it very difficult to link the observed neutrino signal to 
ZAMS mass without additional constraints from classical as- 
tronomical observations in the electromagnetic spectrum. 

In the discussion of SN 1987 A, we have ignored uncertain- 
ties regarding nuclear EOS, rotation of the progenitor core, 
distance, and collective neutrino oscillations. All affect the 
preexplosion neutrino signature in ways that may be degener- 
ate with variations in the compactness parameter. However, as 
we have shown, large variations in the stiffness of the nuclear 
EOS in the density and temperature regime relevant in the pre- 
explosion phase can be disentangled from via observations 
of the emitted neutrino spectrum. Rapid rotation, which de- 
creases both (angle-averaged) luminosities and average neu- 
trino energies, can be constrained by gravitationa l wave ob- 
servations of a galactic event (e.g.,|Ott et al.|2012|. Distance 
uncertainties from electromagnetic observations, which trans- 
late into uncertainties in the absolute luminosities, can be re- 
duced by exploitin g the generic neutr onization burst as a stan- 
dard candle (Kach elrieB et al.|2005| ). Collective neutrino os- 
cillations are not yet fully understood and have not been di- 



rectly incorporated in neutrino radiation-hydrodynamics sim- 
ulations. They may lead to one or multiple energy-dependent 
swaps of spectra between flavors and could thus complicate 
the mapping between neutrino signal and compactness pa- 
rameter. However, the recent understanding suggests that col- 
lective oscillations may not be significant in the preexplo sion 
phase (Chakraborty et al. 2011b a; SarikasJF^|2ofe} but: 
|Cherry et al.|2012| and |Dasgupta et al.|2012|| ). 



Our goal with this study was to highlight overall trends 
of the preexplosion neutrino signal with progenitor structure 
in the limit of spherical symmetry. We expect these overall 
trends to be robust and to carry over to the multi-dimensional 
case. We did not aim at making precise and robust quanti- 
tative predictions for any individual model. These are not 
possible with our current spherically- symmetric radiation- 
hydrodynamics treatment, which neglects inela stic scatter- 
ing, red shift and velocity dependent terms (e.g., |Lentz et al. 
|2012a|b) ). Rotation is included only in an angle- averaged 
manner and other important multi-dimensional dynamics, in 
partic ular convec tion and the standing-accretion shock insta- 
bility (|B uras e t al. 2006; Ott et al. 2008; Ma rek & Janka 
f2009; Marek et al. 2009; Lun d eral.|201 0; Bran dt et al.|2011t 
|Ott et al. 20 12), cannot be captured. Future work to remove 
these limitations will be necessary to produce the reliable sig- 
nal predictions necessary for drawing quantitative conclusions 
from neutrino observations of the next nearby core-collapse 
supernova. However, even with fully realistic and complete 
simulation codes, large parameter studies will be necessary to 
account for prevailing uncertainties in the nuclear EOS and/or 
neutrino interaction physics. 
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